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We study fluctuations of the Wigner time delay for open (scattering) systems which 
exhibit mixed dynamics in the classical limit. It is shown that in the semiclassical limit the 
time delay fluctuations have a distribution that differs markedly from those which describe 
fully chaotic (or strongly disordered) systems: their moments have a power law dependence 
on a semiclassical parameter, with exponents that are rational fractions. These exponents 
are obtained from bifurcating periodic orbits trapped in the system. They are universal in 
situations where sufficiently long orbits contribute. We illustrate the influence of bifurcations 
on the time delay numerically using an open quantum map. 

§1. Introduction 

In closed systems, for example closed billiards, quantum energy levels and eigen- 
functions are related to classical periodic orbits in the semiclassical limit via trace 
formulae. In the case of the energy level spectrum, the trace formula was first devel- 
oped by GutzwiHer. 1 ) Combining trace formulae with statistical information about 
the classical dynamics (e.g. ergodicity) underpins the semiclassical theory of quan- 
tum fluctuation statistics in closed systems. For example, this approach forms the 
basis of attempts to understand universality and the connection with random matrix 
theory in the spectral statistics of fully chaotic systems, and likewise the connection 
with Poisson statistics in regular systems. 

In closed systems that have a mixed (i.e. partly regular and partly chaotic) 
classical limit, some quantum fluctuation statistics are well described semiclassically 
by weighted averages over the regular and chaotic components. However others, 
for example moments of the fluctuations of the spectral counting function around 
the Weyl mean, or moments of the fluctuations of eigenfunctions about the quantum 
ergodic limit, are dominated in the semiclassical limit by classical periodic orbits that 
are close to bifurcation, where their contribution to the trace formulae is enhanced 
by a power of Planck's constant that depends on the nature of the bifurcation in 
question. 2 )~ 5 ) When long orbits contribute it is assumed that one can average over all 
of the generic bifurcations. There is then a competition as to which dominates the 
semiclassical moment asymptotics: essentially the more complicated a bifurcation 
is the larger is its contribution but the smaller is its range of influence. In the 
cases studied so far 3 )' 4 - 1 it happens that for any given moment there is a bifurcation 
which dominates the competition. Different moments are dominated by different 
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bifurcations. The result is that each moment scales semiclassically as a power of 
Planck's constant. These exponents are universal. They take the form of rational 
fractions whose values are given by simple formulae. This differs markedly from the 
behaviour one sees in either fully chaotic or fully regular systems. 

Open (scattering) systems also exhibit universal quantum fluctuation statistics, 
for example in the conductance and the time delay, that depend on whether the 
classical dynamics is regular or chaotic. 6 -* In chaotic scattering systems the fluc- 
tuation statistics coincide with those of random matrix theory. The question then 
arises as to what happens in open systems in which the classical dynamics is mixed: 
is there any analogue of the bifurcation-dominated fluctuation statistics found in 
closed systems? One might initially think not, because the classical trajectories that 
underlie the semiclassical expression for the 5-matrix are scattering orbits, not pe- 
riodic orbits. However, it turns out that some quantities related to scattering can 
be re-expressed directly in terms of the periodic orbits trapped inside the scatter- 
ing region. For example, if the scattering system is a billiard with holes cut in the 
perimeter, these quantities can be expressed either in terms of orbits entering and 
exiting through the holes, or in terms of periodic orbits that never hit the holes. 
Thus in mixed open systems bifurcating periodic orbits can semiclassically dominate 
the related fluctuation statistics in exactly the same way as in closed systems, giving 
rise to new classes of universal scaling exponents that cannot be described within a 
random matrix model. One example where this was recently shown to be the case 
is the conductance fluctuations in antidot lattices. 7 ) 

Our purpose here is to point out that the Wigner time delay is another gen- 
eral example. Specifically, we argue that in mixed open systems fluctuations in the 
time delay may be dominated by classical periodic orbit bifurcations, and that when 
long orbits contribute the fluctuation moments scale semiclassically with universal 
exponents whose values, again rational fractions, are related directly to those calcu- 
lated previously. We illustrate our theory with numerical computations for a class 
of quantum maps. 

§2. The Wigner time delay 

The concept of time delay in quantum scattering was introduced by Eisenbud 8 ) 
and Wigner 9 ) in the context of one-channel spherical wave scattering. Later, Smith 10 ) 
extended the notion to the M-channel case by introducing the lifetime matrix 

M , 

Qij = - ih J2sl—S cj (E), (1) 

c=l 

where S is the standard scattering matrix and the sum runs over all M open channels 
denoted by c. The time delay is defined to be the average of the eigenvalues of Q: 

, . ih i d „ ih d , 

T (E) = -—trS f — S = - — — logdetS". (2) 
M dE M dE v ' 

It can be interpreted as the typical time spent by a scattered particle in the interac- 
tion region. 
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The time delay turns out to be very closely related to the density of states in a 
closed system. This allows one to use the well developed semiclassical apparatus for 
the density of states of closed systems to unravel features of open ones. The con- 
nection was identified (independently 11 )) by Friedel 12 ) and Lifshitz: 13 ) the average 
time delay ([2]) is related to the difference between the level density of an interacting 
Hamiltonian H with respect to a free or reference Hamiltonian Hq. 

Friedel's formalism was used by Balian and Bloch 14 ) to derive a semiclassical 
expression for the time delay which is closely related to the Gutzwiller trace formula 
for the density of states. This splits into a smooth part t(E) and a fluctuating part 
t^{E) determined by periodic orbits: 

T(E) = 2 ^p(E)* 2 -^(p(E)+P & (E)), (3) 

where p is the renormalized density of states. It has a natural interpretation in the 
context of inside-outside duality. 15 )' 16 ) 

The smooth term on the right-hand side of ([3]) can be interpreted as the mean 
density of scattering resonances. It represents the mean time spent in the scattering 
system. The fluctuating term is given by 

p\E) « J- Re £ f; T ~< exp (A- £mA (4) 

* h Tm^i yj\ det(M7-/)| \ h 2 J 

where the sum runs over the periodic orbits 7 which are trapped in the repeller and 
their repetitions. T 7 , <S 7 , M 7 , and p y are respectively period, action, stability matrix 
and Maslov index of the orbit 7. 

The semiclassical analysis of the time delay is often based on the relation to 
the trapped periodic orbits. For example, in fully chaotic open systems, statistical 
properties of these orbits can be invoked to justify the use of random matrix theory 
to model the fluctuations. 16 )~ 18 ) It is worth pointing out that all the emphasis so 
far has been placed on systems in which the classical limit is fully chaotic. This is 
not however the typical situation in actual experiments, where real cavities or soft 
potentials lead to mixed dynamics. Our aim here is to go beyond this idealization by 
incorporating the contributions from orbit bifurcations. This is important, because 
the contributions are, in certain regimes, actually the dominant ones. 



§3. The time delay for quantum maps 



We illustrate the general theory of the last section using quantum maps, since 
in this case the relation between the time delay and the trapped periodic orbits can 
be derived straightforwardly. 

We start with a map that acts on a phase space corresponding to the unit torus. 
Its quantization is defined by a unitary time evolution operator U with finite dimen- 
sion N. The energy dependence of Hamiltonian systems is simulated by including 
a phase factor U = e l£ U that depends on the quasi-energy e. The map is then 
opened up by removing vertical strips from phase space. These play the role of 
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holes in the boundary of open billiards. The corresponding scattering matrix S is 
an M-dimensional matrix, where M < N is the total number of position states in 
the opening. It is the unitary matrix for the transition from the opening onto itself 
after an arbitrary number of iterations of the internal map 19 ) 

S(e) = U 00 + Y1 UoiUhUio ■ (5) 

n=0 

The quantities Uki with k, I G {0,1} are restrictions of the evolution operator to 
the inside and/or outside. They are defined by Uki = PkUP^ ', where Po and Pj are 
the projection matrices onto the opening and the interior, respectively. Po and Pj 
have dimensions M x N and (N — M) x N, respectively, and Pq and Pj are the 
corresponding transpose matrices. The projection operators satisfy PqPo + Pj Pi = 
In, PoPo = ^ M > an d PiPf = In-m where II denotes the L x L unitary matrix. 
The sum over n in © can be performed and the S-matrix written as 

S(e) = P V-—UPI . (6) 

I N - UPjPj 

Planck's constant for quantum torus maps is given by h = (27riV) -1 , and so the time 
delay has the form 

He) = 1 4- log det S. (7) 

In order to derive a formula for the time delay in terms of the periodic orbits of 
the map, it is useful to apply an identity for determinants due to Jacobi. Let A be 
an iV x N matrix given in terms of the auxiliary block matrices B,C,D and E, and 
A^ 1 its inverse matrix given in terms of W, X, Y and Z: 

B and W are assumed to be square matrices with the same dimension. Jacobi's 
determinant identity then states that det B = det Z det A. 

Let us identify A with the N x N matrix [Jjy — IJ Pj P{\~ X IJ and the subblock 
B with PoAPq = S. Then Z follows as PjA~ x Pj and the Jacobi identity leads to 

det S = det([Jjv - UPfP^U) det(P/^ 1 [/ 7V - UPjPi]Pj). (9) 
A little elementary linear algebra then gives 



detS = detU detiUlI - lN Z M) . (10) 
det(/Ar_ A / - Un) 

For the evaluation of the derivative of the logarithm of det S we need 

— log det [7= — log e i7Ve det [/ = iiV , (11) 
de de 
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and 



d_ 

ds 



oo 



logdet(ijv-Af - Uii) = 



e m£ tr U\ 



■n 



(12) 



n=l 



This leads to the final result 




MN \ 2vr 



1 / N 




e ine tiU\ 



n 



) 



(13) 



which is in agreement with the general formula ([3]). The first term on the r.h.s. of 
(|13|) represents the mean time spent in the scattering region, (r(e)), while the second 
term is the fluctuating part of the time delay T fl (e). It should be noted that r fl (e) 
is written in terms of Uu rather than U, the evolution operator for the unopened 
map. The powers of traces of Uji are semiclassically related to the periodic orbits 
which lie completely in the interior of the open map. If the traces are evaluated in 
the semiclassical approximation one reproduces ([3]) and (JH). 



Generic systems have a phase space in which regular islands and regions of 
chaotic sea coexist. One of the main characteristics of mixed systems is the bi- 
furcation of periodic orbits. Bifurcations are events where different periodic orbits 
coalesce when parameters of the system are varied. They are important in semiclas- 
sical approximations because bifurcating orbits carry a semiclassical weight that is 
higher than that of the isolated (unstable) periodic orbits and sometimes even that 
of tori of regular orbits. 

Consider the Gutzwiller contribution of a periodic orbit 7 and its repetitions m. 
The amplitude in diverges if has an eigenvalue one, which happens at a bi- 
furcation. This is because periodic orbits are assumed to be isolated in the derivation 
of the trace formula. 1 ) More specifically, the trace formula can be derived by inte- 
grating over Poincare sections perpendicular to periodic orbits, and periodic orbits 
appear as stationary points of these integrals. Consequently, the usual stationary 
phase approximation breaks down when stationary point coalesce. The remedy is 
to perform a uniform asymptotic expansion valid throughout the bifurcation pro- 
cess. 20 ^ 24 ) This is obtained by rederiving the trace formula using the appropriate 
generating function for the Poincare map <&{Q , P) from (Q, P) to {Q , P ) in normal 
form coordinates. For two-dimensional systems the semiclassical contribution of a 
bifurcation to the density of states is then given by 



As an example we take a saddle node bifurcation, <P(Q, P) = P 2 + x\Q + Q 3 , where 
x\ depends on the energy or other system parameters and vanishes at the bifurca- 
tion. The stationary points occur for negative x\ at (Q,P) = (±y/—xi, 0). The 
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contribution (|14h at the bifurcation is oc /i _/3 , where j3 = 7/6. Away from the bi- 
furcation (i.e. when —xi/H is large) one obtains contributions of isolated orbits that 
are oc h" 1 . The two regimes are interpolated by an Airy function. 

Besides (3 there are further exponents that are important for the semiclassical 
influence of the bifurcation. They describe the size of the parameter intervals over 
which the bifurcation is semiclassically stronger than isolated periodic orbits. Con- 
sider <P(Q, P) = P 2 + xiQ + Q 3 in the example. We can make the exponent in the 
integral (|14p ^-independent by scaling Q = Qh 1 ^ 3 , P = Ph 1 / 2 and x\ = x\ 
Hence the relevant x\ interval scales like h ai , where o\ = 2/3 in this example. The 
same analysis can be applied to more complicated bifurcations which have more 
parameters Xj, i = 1, . . . , K, in their normal form. (K denotes the codimension of 
the bifurcation.) One then obtains a characteristic exponent Oi for every parame- 
ter Xj. Because of this finite extension in parameter space, bifurcations of higher 
codimension have to be taken into account even if only one parameter is varied. 

Generic bifurcations of periodic orbits are characterized by the codimension K 
and the repetition number m of the orbit for which the bifurcation occurs. A sys- 
tematic investigation of the influence of the different bifurcations on moments of 
the density of states was carried out by Berry, Keating and Schomerus. 3 ) This has 
subsequently been extended to determine their influence on the statistics of wave- 
functions 4 )' 25 - 1 and, more recently, on the moments of the conductance fluctuations 
in antidot lattices. 7 ) 

Fluctuations of the Wigner time delay can be characterized by their moments, 
defined as 



where {•••)e,x denotes averaging over energy and over parameter space. In any 
parameter interval of a system with mixed dynamics infinitely many bifurcations 
occur, most of them for very long periodic orbits. If h is small enough, then these 
bifurcations are important. To determine the influence of a particular bifurcation 
with codimension K and repetition number m on the moments of the time delay, 
we replace the average in (|15p by an average over parameters in the normal form. 
Performing the scaling procedure 3 ) we can extract the h dependence for the different 
bifurcations M 2 k,m,K ~ K~ nk < m < K where r]k, m ,K = 2k(3 myK - J2iLi &i,m,K- Hence we 
see that the importance of a bifurcation depends on the quantity r\k,m,Ki which in 
turn depends on the characteristic exponents f3 m ,K and crj, mj x- The bifurcation that 
is most important for a particular moment is that for which rjk,mK is largest, and 
hence one finds that A^fc ~ /t _r?fe where % = max my K (Vk,m,K)- There is a different 
winner of this competition between bifurcations for every k. It follows from ([2D 
that all the different exponents for the time delay (|15p coincide with those for the 
corresponding moments of the density of states. 3 )' 7 - 1 




(15) 



Periodic orbit bifurcations and scattering time delay fluctuations 7 

§5. Numerical results 

We now illustrate some of the general ideas described above using a family of 
perturbed cat maps. 26 )' 27 ) These are maps of the form 

CMs OCH^GH 1 (16) 

where q and p are coordinates on the unit two-torus, and are taken to be a position 
and its conjugate momentum. We will concentrate on one particular bifurcation and 
investigate its influence on the second moment of the time delay. Using (|15|) , ([3D and 
(p~3|) one can write 

M 2 =± f\p fl A {e)?de^^Y,\^\ 2 e- n2A \ (17) 

where p\ is the fluctuating part of p convoluted with a normalized gaussian of width 
A. In our computations A was taken large enough so that the dominant contributions 
to p* come from the n = 1 term in the sum in (|17p . 

For k = the map (|16p is uniformly hyperbolic. The perturbed map is guaran- 
teed by Anosov's theorem to be strongly chaotic for k < (y/3 — l)/v5 ~ 0.333. The 
period-1 fixed points at 

Qj = 2 ' ~ ~h cos ( 27rqj ^) ^ 

for integers j=0,l are then unstable. 27 ) Outside this parameter range bifurcations 
can occur leading to a mixed phase space. When k = k^h = 5.94338 the phase space 
is almost entirely ergodic but a saddle node bifurcation gives rise to a new pair of 
period-1 orbits. 2 ) For this map, it has been shown 27 )' 28 ) that trU can be expressed 
in the form 

tvU= y/(N/i) / exp^niNSjiq^dq (19) 

where Sj(q) = q 2 + sin (2i:q) — jq, so that the phase is stationary at the fixed 
points (PTBj) . 

After removing strips from phase space one gets the open map. In all the cases 
illustrated in this section, the ratio M/N of the dimension of the quantum map to the 
number of open channels is 0.28. In order to see the contribution of bifurcating points 
to the moments (|17p we have computed |trt///| 2 in two different configurations. 
First, two strips are located in such a way that they block the unstable fixed points 
go = 0.81425 and q\ = 0.6857. Second, the two strips are placed such that they block 
the bifurcating points q^ = 0.4453 and qf = 0.0548. These are period-1 fixed points 
which originate in a saddle node bifurcation at Kbif- The logarithm of the second 
moment A^2 is plotted in figure[T]as a function of the perturbation parameter k. The 
two strips block either the bifurcating points (solid black curve and crosses) or the 
isolated points (dashed red curve and circles). In this and the subsequent two figures 
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Fig. 1. log M2 versus k with A = 1.5 and N = 997. Solid black curve and crosses for strips centred 
at the bifurcating fixed points; dashed red curve and red circles for strips centred at the isolated 
fixed points. Curves and symbols for theory and numerics respectively. 



the theoretical curves correspond to asymptotic evaluations of (fTUj) . 2 ) as described 
in outline in the previous section. 

Figure [2] illustrates the contribution of isolated points to M2 as a function of 
the dimension N. In this case. striDS are centred at the bifurcatine Doints so that 
only the isolatec ^ 




ln(N) 



Fig. 2. log M2 versus log N for A = 1.5 and k = ftbif. Strips centred at the bifurcating fixed points. 
Solid black curve and crosses for theory and direct numerical evaluation of the trace using (|17[l 
respectively. 

By contrast, figure [3] illustrates the situation where the strips block the isolated 
orbits. The contribution of the remaining, bifurcating orbits to | tr U\ 2 grows semi- 
classically like iV 1 / 3 in this case 2 ) (i.e. the scaling exponent is 1/3). For k around 
Kbif, as in figure [H tvUji is more intricate and involves Airy functions. 2 ) 

These figures illustrate clearly how the bifurcating trapped periodic orbit dom- 
inates the fluctuations of the time delay semiclassically. 
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. 3. \0gM2 versus logiV for A = 1.5 and k — Kbit- Strips centred at the isolated fixed points. 
Dashed red curve and circles for theory and direct numerical evaluation of the trace using (|17|l 
respectively. 
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